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Within the framework of two-flavor covariant baryon chiral perturbation theory we have expressed 
the Chernyak-Zhitnitsky, Ioffe and Dosch currents in terms of chiral fields to provide leading one- 
loop extrapolation formulae for the leading and next-to- leading twist normalization constants /jv, 
Ai and A2. Finite volume effects due to pion loops have been taken into account. The occurring 
low energy constants are fitted to data obtained from recent lattice QCD simulations in order to 
extract the values at the physical point. 

I. INTRODUCTION 

In recent years, the combination of lattice quantum chromodynamics ((L)QCD) and effective field theories (EFTs) 
has proven to be especially fruitful. EFTs allow to relate hadronic properties in the real world to those in a fictitious 
one, where the quarks are heavier than in reality and are restricted to a finite volume (at least for not too small 
volumes). Hadron properties in such a theory are much cheaper to calculate numerically and this is what is usually 
done. This work addresses one specific item in this general broad program, namely the quark mass and volume 
dependence of the zeroth moments of specific distribution amplitudes (DAs) , which are needed for the description of 
exclusive processes. 

The nucleon DAs are non-perturbative functions describing the internal nucleon structure. They are most conveniently 
defined by nucleon-to- vacuum matrix elements of nonlocal light-ray three-quark operators |lH3( . Restricting ourselves 
to contributions of (next-to-)leading twist and leading conformal spin, only local operators without derivatives have 
to be taken into account. Therefore only three independent structures remain, involving three (scale-dependent) non- 
perturbative constants /n, Ai and A2, called leading and next-to-leading twist wave function normalization constants. 
Information on nucleon DAs can not be obtained in a direct way from experiments. They enter observables only 
through convolution integrals and theoretical guidance is needed to come up with sensible parametrizations. There 
are two theoretical frameworks capable of extracting information on DAs, which are of completely non-perturbative 
nature: the analytic QCD sum rule approach and the numeric approach of LQCD (see e.g. Ref. Q). The QCD 
sum rule technique is well established and leads to deep qualitative insights. However, it suffers from big systematic 
uncertainties, typically of order 10% — 30%. LQCD can provide high precision results, but only if all systematic 
uncertainties are well under control, including the quark mass dependence and finite volume effects. 
The goal of this work is to set the chiral extrapolation of the wave function normalization constants mentioned above 
on a firm basis by the use of covariant baryon chiral perturbation theory (BChPT). The construction of the nucleon 
interpolating currents that project onto the normalization constants in terms of chiral fields forms the basis of this 
analysis. The currents corresponding to the next-to-leading twist normalization constants Ai and A2 appear also 
in grand unified theories, where they mediate proton decay, and have been constructed for the three-flavor case at 
leading order by Claudson, Hall and Wise in Ref. [Sj, which served us as starting point. Using the same technique we 
construct the three interpolating currents corresponding to the wave function normalization constants in a two-flavor 
framework and calculate the respective matrix elements to leading one-loop order, where we take into account finite 
volume effects due to pion loops. 

This work is organized as follows: in Section |H] we give a brief introduction to chiral perturbation theory (ChPT) and 
to nucleon interpolating currents. The construction of the currents in terms of chiral fields is performed in Section Hill 
In Section IIVI we evaluate the contributing loop diagrams to obtain extrapolation formulae for the wave function 
normalization constants and fit the occurring low energy constants (LECs) to recent QCDSF lattice data in order to 
determine the values for physical masses. We conclude and give an outlook in Section IVl 

II. THEORETICAL SETUP 
A. Covariant baryon chiral perturbation theory 

There exist quite a number of formulations of effective hadronic theories. In this work we use covariant BChPT 
for reasons of convenience: manifest Lorentz covariance simplifies the construction of the low energy form of nucleon 
interpolating currents and allows a direct comparison to the general decomposition of the proton-to- vacuum matrix 
elements given in Ref. [2^ . 
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Chiral effective field theory (ChEFT) is based on the spontaneous breakdown of chiral symmetry in QCD. The 
additional (small) explicit symmetry breaking induced by finite quark masses yields massive quasi Goldstone bosons 
which can be interpreted as mesons. By establishing a connection between the set of left cosets of the unbroken 
vectorial subgroup in the chiral group (G = Sf7(n/)i <8> SU{nf)n) and the Goldstone boson field configurations one 
finds that the mesonic degrees of freedom can be represented by a SU(rif) matrix U which transforms linearly under 
a chiral rotation g — (L,R) £ G (see Refs. @, 0|) : 

U ->■ RUL^ , (1) 

i.e. the behavior of the Goldstone boson fields under symmetry transformations is determined by geometry (up to 
reparametrizations). It is most convenient and common practice to choose canonical coordinates for the parametriza- 
tion of U [3l|. In the case of two light quark flavors (n/ = 2) U then reads: 

U = e^* with $ ee ttV 1 , (2) 

where is the pion decay constant in the chiral limit, a 1 (i — 1,2,3) are the Pauli matrices and it 1 (i — 1,2,3) 
represent the pion fields. The nucleon fields (proton and neutron) are grouped in a so-called isospinor 

*-f;V (3) 



Its transformation properties under the chiral group are only fixed on the vectorial subgroup where it has to transform 
as a doublet, i.e. we are free to choose the most convenient realization of chiral symmetry from various (physically 
equivalent 0, H[) possibilities. Therefore we define the field 



u ee VU = * , (4) 

which has to transform under chiral rotations as 

u -> RuK^ = KuL f (5) 

in order to reproduce the transformation properties of U . Eq. ((3J) defines the compensator field K. Since K meets 
the condition stated above, we can choose the nucleon field to transform as 

$ -> . (6) 

This choice allows for a convenient construction of invariants because one can write down mesonic building blocks 
that transform as X -> KXK\ namely |32||: 



i [u 



X± 



= .,t 



tfyu - ufyut) , (7) 
v)x^ ± ux^u , (8) 



where x = ^B^M includes the quark mass matrix M = diag (m u , m^) and the condensate parameter Bo- The 
price one has to pay for adjusting the transformation properties of baryonic and mesonic fields this way is that 
non-vectorial chiral transformations are local since they depend on the mesonic field u. Hence one has to introduce 
covariant derivatives: 

ee (0„ + r M ) * , (9) 

D^X = d^X+[T^X] , (10) 
where T M is called the chiral connection and is given as 

r„ = i (vld^u + ud^) . (n) 

In Table U we have summarized the building blocks together with their main properties for convenient reference. In 
the baryonic sector the nucleon three-momentum and the pion four-momentum (including the pion mass) serve as 
(small) parameters of the perturbative expansion generically denoted as p. The power of p in a term defines its chiral 
dimension D. The lowest order pion- nucleon effective Lagrangian reads [9( 



4£r =1) = * (0 - m + ^75) * , (12) 

where is the axial coupling constant and rop is the nucleon mass (both in the chiral limit). For a full list of terms 
up to fourth order we refer the reader to Ref. [10( . To allow for both a systematic expansion of loop contributions in 
powers of p and manifest Lorentz covariance we apply the infrared regularization (IR) scheme [ll| . 
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TABLE I: Chiral dimension and transformation properties of the basic fields. Note that the chiral dimension of a derivative 
depends on the field it acts on: it has dimension zero/one when acting on a nucleon/pion field. 







u 


tV 


X± 


Chiral dimension 








1 


2 


SU(2) L ®SU(2) R 




RuK 1 = Kutf 


Ku^ 


K X± K^ 


Parity 


70* 


«t 




±X± 



B. Three-quark interpolating currents 



We are interested in proton-to- vacuum matrix elements of three-quark currents. A general decomposition for such 
overlap matrix elements was given in Ref. Restricting ourselves to contributions of (next-to-)leading twist and 
leading conformal spin, we can reduce the problem to the case where all quark fields are located at the origin. For 
this situation one obtains a simplified (local) decomposition of the matrix element: 

(0\e^ k uiu^\N(k)) = Vf($C) a0 ( l5 N) 7 + Vlm^CU^N^ 

where \N(k)) represents a proton state with four-momentum k (k 2 = m 2 N ), N denotes the corresponding spinor and 
u, d are quark field operators. The latin letters refer to color and the greek letters to Dirac indices. C denotes the 
charge conjugation matrix. This decomposition has three independent parameters, since = 7~i due to isospin 
constraints [2]. In order to project onto the corresponding components we introduce the three standard baryonic 
currents: the Ioffe (I), Dosch (D) and Chernyak-Zhitnitsky (CZ) current. They have been defined in [l2T4l4j as 



V \x) = e^ k (u <r (a:)C77 M «»'(i)) 757^ fc (x) , (14a) 
V D {x) = e^ k (u lT (x)Ca^uHx)^ l5 a^d k (x) , (14b) 
V cz (x) = ^e ijk [{u lT (x)C^(x)) l5 ^d k (x) - (u tT (x)Cjtd^(x)) ^u k (x)] , (14c) 



where z is an arbitrary, dimensionless, light-like auxiliary vector. Note that the given definition of the Chernyak- 
Zhitnitsky current is not the original one but an isospin improved version which does not contain A + baryon pollutions. 
The prefactor has been chosen such that the original current normalization is maintained. Using the decomposi- 
tion (|T3]) one obtains: 

{0\ri l {Q)\N{k)) = (V? - 4V^) m N N = X imN N , (15a) 
<0|77 D (0)|iV(fc)) = 6 (Ti - 4T 3 °) m N N = \ 2 m N N , (15b) 

(0|r, cz (0)|iV(fc)) = 1 (V? + 2T?) (z ■ k)jN = f N (z ■ k)jN . (15c) 



III. CONSTRUCTION OF CURRENTS 



In the following we will construct the currents given in Eq. (fl"4f in terms of the building blocks shown in Table Q] 
For that purpose we analyze the symmetry properties of the currents under chiral rotations, parity transformation 
and (for the Chernyak-Zhitnitsky current) the separation in plus and minus components. At least for the two-flavor 
case the Ioffe and Dosch currents exhibit similar behavior under chiral transformations and thus can be treated within 
one approach. We refer to [l5j for a comprehensive presentation of the chiral transformation properties of three-quark 
operators. 
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A. Ioffe and Dosch current 



Using the method of Fierz transformation [l|| [TtJ one finds that the currents can be separated in parts of different 
helicity. Defining q p = (1, 0) T to project onto the quark content of the proton we can rewrite the currents as 

T 1 \x)=^( lR U L {x)- lL U R {x)) , (16) 
Tf{x) = %ql( lR U R {x)- lL U L {x)) , (17) 



where we have defined 

u l,r [x) _ jjk ^ T L/R (x)C(iT 2 )q^ L/R (x)) q\x) . (18) 

q = (u,d) T is a quark doublet field and qi/u = Jl/r ? = (1 T 7s)/2 q its left- and right-handed projection, r 2 
is the second generator of the SU(2) Lie algebra in the fundamental representation. By applying an infinitesimal 
chiral transformation to the operators ^lU L ^ R and ^rU L ^ R one finds that they transform as (2^,1^) and (1^,2^) 
respectively. Hence, we can write the low energy form of the currents as 



71=0 fe=l 
OO i n 

EEa 



LR ' 



(ra) / n (n) _ n (n) 
k \ w k,RR ^kXL) I 



n=0 fc=l 



(19) 
(20) 



where we have suppressed the x dependence, n gives the chiral dimension of the operator and i n is the number of 
operators with dimension n, which are labeled by k. Due to their similar transformation properties the Ioffe and 
Dosch currents only differ in the prefactor and the LECs , (3^ , i.e. 0^ RL = 0^' RR and 0^' LR — 0^ LL . The 
different current operators that reproduce the transformation properties given above are listed in Table HH sorted by 
their chiral dimension. The different sign in some of the terms ensures correct (positive) parity. The terms which 
contribute to the proton-to-vacuum matrix element at leading one-loop order (which equals second chiral order) are 
marked by (*). 



TABLE II: Low energy operators for the Ioffe and Dosch current grouped according to their chiral dimension n. The list of 
second order terms is far from complete but solely contains terms that contribute to the proton-to- vacuum matrix elements of 
the currents at leading one-loop level (indicated by (*)). The contributions which vanish in the limit of exact isospin symmetry 
are marked with (7). 



o (n) - o K 

k,RL — k,RR 



-,(») 



o (n) - o y 

k,LR — k,LL 



wyR^! 

u (X+ - ^tr{x+»7fl* 



— u'jt^ 

" t M, l 7L7 M * 
u t M M 7Li-D M * 

-tr{x+}«V* 
V (x+ - I tr {x+}) 7i* 



(*) 



(*) 



B. Chernyak-Zhitnitsky current 



For the construction of the current the transformation properties under chiral rotations play a key role. Splitting 
up the current into parts of different helicity by the use of Fierz transformation one obtains 

r 1 CZ {x) = ^T a ( 1 dA R > a {x)- lR U L ' a {x)) , (21) 
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where we have defined 

Ql,r^ {x) = ^ T L/R ( x )Ca^(tT 2 )T a q\ /R (x)) a;q k (x) . (22) 

Applying an infinitesimal chiral transformation one finds that the diquark part does not transform in the trivial 
representation (as it was the case for the Ioffe and Dosch current) but under the three-dimensional adjoint represen- 
tation. Therefore the operators jL^ R ' a and r ynU L a transform as (2l,3r) and (3l,2r) respectively. Suppressing the 
x dependence the low energy version of the Chernyak-Zhitnitsky current takes the form 

v cz - E E 4 n) (0£is - o£&) . (23) 

n=0fc=l 

The operators oI^'r and O^m, are listed in Table IIIII The relative sign is fixed by the condition that the total 
Chernyak-Zhitnitsky current has positive parity. The absence of terms with two derivatives acting on the nucleon 
field ((z ■ D) 2 ^) needs some further comment. Consider the plus/minus projection operators which have been defined 
to separate components of lower and higher twist of a nucleon spinor with momentum k (see e.g. Ref. Q): 

which, in position space, read 

A + ,^, A_ S ^g. (25, 

2z ■ a 2z ■ a 



Using the definition given in Eq. (|14c|) it is easy to verify the following properties of the Chernyak-Zhitnitsky current: 

A +V cz = , A_ ?/ cz = rf z . (26) 

For an application of the plus/minus projectors on the low energy version of the current recall that the chiral connection 
is of first chiral order, since it contains a derivative acting on a meson field. Accordingly, we can replace the partial 
derivative by a covariant one (the difference is of higher order). Hence, we find that the minus projection of terms 
including two covariant derivatives acting on the nucleon has the form jt(z ■ D)^f: 

A_(z • D)H = ^(z- D) 2 ^ = -p-{z ■ D) 2 * + 0(p) ex fa ■ D)* + 0(p) , (27) 
2z ■ a 2z ■ D 

what justifies neglecting these terms in the operator. It is worth noting that for terms proportional to £{z ■ D)^ we 
have the exact relation 



K^{z ■ £>)* = ^t(z D )*=[ 1 -2^d) ^ z ■ = & z ■ D )* ■ ( 28 ) 
A similar argument leads to the structures jtz^ and f{z ■ D)^^ for the terms of first chiral order. 

IV. RESULTS 

A. Extrapolation formulae 

Dimensional counting yields that one has to consider the Feynman diagrams shown in Fig. [T] to calculate the matrix 
elements given in Eq. (fT5j) at leading one-loop level. The squares depicting the operators have to be derived from 
Eqs. (flT))) . (|2T)|) and (|2"3")l for the respective current, while the circles denote standard vertices originating from the 
effective Lagrangian. The digits within the vertices and operators give their chiral dimension. The leading tree- 
level diagrams come with a Z-factor contribution. Since we work in the limit of exact isospin symmetry, the terms 
marked with (J) in Tables [TTI and HITl do not contribute. The results for particular diagrams can be found in App. [Cl 



Interpreting the contribution of Fig. 1(b) as a chiral counterterm, i.e. splitting the corresponding LEC into a finite 
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TABLE III: Low energy operators for the Chernyak-Zhitnitsky current grouped according to their chiral dimension n. The list 
of second order terms is far from complete but solely contains terms that contribute to the proton-to-vacuum matrix element of 
the currents at leading one-loop level (indicated by (*)). The contributions which vanish in the limit of exact isospin symmetry 
are marked with (7). 



n 


k 


s-\{ n) ,a 

k,LR 


s~\(n) ,a 
°k,RL 







1 


z M z v u 1 " 2 ^u-jL-f^iD^^ 


-z ,i z l, u 2 T a u t 'y R "f l _ l iD u i i! 


(*) 


1 


1 




2* 'Z 1 'uUp,UT a U ^jR'fv^ 




1 


2 


Z Z U U U T U'yL^n'yulDvW 


Z^Z UU UT u'^Rj^al-Uv W 




1 


3 


Z^ Z V V^ 2 T a UU^L^jv^ 






1 


4 


z tl z u u x2 T a uu a ^ L ^^ a iD v ^ 


2 M 2 v u 2 T a u* ti Q 7fl7 fl 7c i£>„ * 




2 


1 


z^z v tr {x+} u t2 r a W7i7 ft i-D I/ * 


-z^z" tr {x+} Ji 2 r a u t 7 fl 7 M iCvE' 


(*) 


2 


2 


(x+ - 5 tr {x+» uVuT^iA,* 


-z^ii (x+ - |tr{x+}) ur a u t 7 ii 7 fl iD y * 


(*),« 


2 


3 


z^u+Vu ( X + -|tr{x+»7i7^* 


-^"uW (x+ - § tr 7k7m^* 


(*),(/) 



(a) 



(b) 



/ \ 
/ \ 
\ / 

\ , , / 







(c) 



k-q 



(d) 

FIG. 1: Feynman graphs relevant for the calculation of current proton-to-vacuum matrix elements to leading one-loop order. 
While (a) is a diagram of zeroth chiral order, (b) (c) and (d) are graphs of second order. 



(denoted by ^) and a divergent part in such a way that the divergence of the loop graphs is cancelled exactly, we 
find the following extrapolation formulae: 



(Aitojv) (m w ) = af' 
(X 2 m N ) (m n ) = 0[ o) ^8 

/^(TOtt) = <4 0) ( 1 



2(4ttF t ) : 

™2 



6 5 i + (3 + 9 ff i)ln. 



(47rf T ) 



8(4^)2 



69 2 a + (3 + 9g\) In ^j) + 32/^'^ + 0(m» ) 
+ (19 + 9g\) In 4 ) ) + 4«i 2) ' (, V + {m*) 



(29a) 
(29b) 
(29c) 



As expected we find an identical pion mass dependence (up to LECs) for Ai and A2 which is not surprising in view 
of the similar low energy form of the Ioffe and Dosch current (cf. Eqs. (flT)|) and ([2H|) L The extrapolation formula for 
/at has a slightly different form as a consequence of a different weighting of the various graphs due to the particular 
form of the current (see App. [C]). Note that the formulae given above are only valid for an infinite volume. For a 
consideration of finite volume effects it is justified to apply the EFT formalism introduced in Refs. [l8l - l2"lj |. since they 
are predominantly caused by mesons traveling around the box (contributions of single quarks walking around the box 
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TABLE IV: QCD sum rule predictions for f N , Xi 
(rescaled using a, one- loop approximation). 


(NLO) and A 2 (LO) at scales of 1 GeV (taken from Refs. @, HI 


) and 2 GeV 




1 GeV 


2 GeV 


fs x 10 3 [GeV 2 ] 


5.1 ±0.8 


5.0 ±0.8 


-A a x 10 3 [GeV 2 ] 


33 ±8 


35 ±8 


A 2 x 10 3 [GeV 2 ] 


51 ± 19 


52 ± 19 



are exponentially suppressed [22]). Therefore we rewrite X e {/at, Aimjv, A2TOat} as follows 

X(m Vl Li = co) = X(m v , L 2 ) + X(m n , Li = go) - X(m„,Li) , (30) 

S ' 

=8X(m^.Li) 

where the part with infinite boxlength is given by the corresponding formula of Eq. (|29[) , while SX can be determined 
by the calculational method presented in Ref. [20| (see App.[B]for the result). Note that by applying this procedure 
we consider the effect of pion loops, but neglect a volume dependence of the effective Lagrangian itself that cannot 
be excluded in our case: in Ref. |2j it was reasoned that for standard boundary conditions (i.e. periodic for bosonic 
and antiperiodic for fermionic fields in space and time direction) the finite and infinite volume effective Lagrangian 
densities resemble each other. However, this condition is not given for the lattice simulation we rely on. Like the 
majority of contemporary LQCD calculations it uses mixed boundary conditions for quark fields (antiperiodic in time 
and periodic in spatial directions). Therefore, strictly speaking, one cannot rely here on the proof given in Ref. (2lT | 
and one cannot exclude that the effective Lagrangian has to be modified by volume-dependent terms. 



B. Fitting of lattice data 

In the following we will determine the occurring LECs by a fit to recent lattice data taken from Ref. 0, [Hj]. 
Two flavors of clover fermions have been used for this lattice simulation with antiperiodic boundary conditions in 
the time direction and periodic boundary conditions in the spatial directions. The values are given at a scale of 
2 GeV and the quarkonium parameter ro = 0.5 fm has been used to determine the lattice spacing. For details on the 
lattice evaluation we refer to Q. In future fits, using new yet unpublished data, QCDSF will use the derived ChPT 
expressions. The comparison presented here does not include these data and is meant as illustration for the general 
trend. Working in a ChEFT framework we have to restrict ourselves to pion masses that can be considered small 
with respect to a typical hadronic scale of 1 GeV. For the given data set the particular value of 600 MeV constitutes 



a justifiable tradeoff between the accuracy of the fit and the reliability of (one-loop) ChPT (cf. Ref. [26j, where an 
upper limit of applicability of BChPT has been estimated to 600 MeV). Table \TV\ contains the QCD sum rule values 
for /at, Ai and A2 at scales of 1 GeV (taken from Refs. @, [23j]) and 2 GeV (rescaled in the one- loop approximation 
with A^g = 226 MeV following Ref. [4] for better comparability). 

Since the extrapolation formulae (|30[) describe the pion mass dependence of the product \i/2m,N one mig ht be tempted 
to divide them by mjv and insert a chiral expansion for the nucleon mass (see for instance Ref. [27|). However, a 
direct fit of the product Xinfnjsf turns out to be superior with respect to both effort and accuracy. The curves in 
Fig. [Dare obtained by fitting the LECs that appear in the extrapolation formulae ([3"U| to lattice data. Since the LECs 
only occur as prefactors, one can use a linear fit procedure. We have chosen the axial cou plin g gA = 1.21, the pion 
decay constant F„ = 86.2 MeV and the nucleon mass mo = 0.88 GeV in the chiral limit [27H29| . The fitted values for 
the relevant LECs are given in Table [V] while extrapolated values of the wave function normalization constants can 
be found in Table IVT1 



V. CONCLUSION AND OUTLOOK 



Within the framework of two-flavor covariant BChPT we have constructed the Ioffe, Dosch and Chernyak-Zhitnitsky 
current. While such a construction has already been performed by Claudson, Hall and Wise for the Ioffe and Dosch 
current in the three-flavor case [5], it is, to our knowledge, entirely new for the Chernyak-Zhitnitsky current. 
By calculating the zeroth moments of the respective nucleon DAs at leading one-loop order we provide chiral ex- 
trapolation formulae for the leading and next-to-leading twist wave function normalization constants /at, Ai and A2. 
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FIG. 2: Fitted curves at infinite volume. The red points represent the original data points, w hile t he blue ones ar e fini te volume 
corrected. The black points at the physical pion mass are the sum rule predictions (NLO for (a) and (b) LO for (c) I and were 
not used for the fit. The colored bands represent the confidence bands for a confidence level of 68.3% and 95.4%. 
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TABLE V: Fitted values for the LECs at leading one-loop order in the limit of exact isospin symmetry and fj. — 2 GeV. 

LEC fitted value 

4 0) x 10 3 [GeV 2 ] 2.5 ±0.1 

nf ] x 10 3 0.6 ± 0.5 

a{ 0) x 10 3 [GeV 3 ] -7.9 ±0.6 

a< 2) xl0 3 [GeV] -10.2 ±1.8 

/3[ 0) x 10 3 [GeV 3 ] 8.0 ±0.6 

f3{ 2) x 10 3 [GeV] 9.8 ±1.8 



TABLE VI: Values for the wave function normalization constants at the physical point. Finite volume effects due to the pion 
cloud have been accounted for. The errors given are taken from the 95% confidence bands. Note that we have neglected the 
systematic error introduced by the quarkonium parameter tq. 



normalization constant 


value at the physical point 


f N x 10 3 [GeV 2 ] 


3.3 ±0.5 


-Ai x 10 3 [GeV 2 ] 


43 ± 10 


A 2 x 10 3 [GeV 2 ] 


86 ± 19 



We took into account finite volume effects due to pion loops. Fitting the occurring LECs to recent lattice data we 
obtained the values at the physical point. The values obtained agree with QCD sum rule expectations except for /jy 
which comes out somewhat small. 

To obtain information on the DAs at non-zero quark separation one has to consider higher moments of the DAs. These 
can be related to matrix elements of local operators containing derivatives acting on quark fields (see Ref. which 
can be evaluated on the lattice. Therefore it is necessary to generalize our results to the case of operators containing 
derivatives in order to provide a systematic chiral extrapolation and finite volume corrections for the respective matrix 
elements. A possible application of the currents constructed in Sec. IIII AI would be the evaluation of nucleon-to-pion 
transition DAs at leading one-loop level. Such a calculation seems promising because the result would depend on the 
same LECs (among others) and the matrix elements are relevant for the main proton decay channel p — > 7r° ± e + in 
grand unified theories. 
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Appendix A: Infrared integrals 

In order to apply the infrared regularization scheme (introduced in Ref. [ll|) one defines the infrared integrals as 
follows: 

7 M 4 - rf f 1 (A1) 

' i J m {2n) d ( m 2 _ q 2 _ le f (M 2 _ (p _ q y _ ie y ' 
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where /i denotes the scale and IR indicates that the Feynman parametrization connecting the meson and baryon 
propagators is extended to infinity. Using the definition 



-1 



(4tt) 2 



1 



1 + ln4?r - 7 B ) 



1 ? 

m m 
■ In — =- 



where e = 4 — d and is the Euler-Mascheroni constant, the solutions to i^o, Io,i and Ii.i read: 

h,o = 2m 2 L 
Io,i = , 

h,i = — —l 



(4tt) 2 



m 2 - M 2 + p 2 



1 



P 



(47T) 



M 2 +p 2 



2p 2 



, m 
l-ln-y 



m 2 (m 2 -Af 2 +p 2 r 
21 ' — ^2 arccos 



- (to 2 - M 2 + p 2 ) 



For an evaluation of the tensorial integrals occurring in our calculations we define 



V 1 k,l ~ 



nv r(2) , jj !/r(3) 



jti 



4-d 
4-d 



IR 



IR 



(27T) d ( m 2 _ q 2 _ le f ( M 2 _ (p _ g )2 _ ie ) 

W (m 2 



q 2 - ief (M 2 - (p - g ) 2 - ze) 1 



(A2) 

(A3a) 
(A3b) 

(A3c) 



(A4a) 
(A4b) 



Identities derived for standard loop integrals also hold for the infrared part, since the decomposition of standard loop 
integrals into an infrared singular part and a regular part is unique. Therefore Lorentz decomposition yields 



(i) 
a-,; 

r(3) 
l k,l 



2^ 



h-i,i + (m 2 - M 2 + p 2 ) I k A 



2p 2 (d 



((2 - d)/ fc -i,i - 2m 2 I k j +d(m 2 - M 2 + p 2 ) l£ 



(A5a) 
(A5b) 



Appendix B: Finite volume corrections 

According to Ref. 20] we obtain the finite volume correction 5X in Eq. (|3"0")) by calculating the difference between 
integration and discrete sum over loop momenta (we neglect a possible volume dependence of the effective Lagrangian 
as mentioned in Sec. IIV A| . Hence, we define 



XT 1 f dq Q J f d 3 q 



2vr [J (2tt) 3 JJi U q J (m 2 - q 2 - ief (M 2 - (p - q) 2 - ief ' ^ 



where Li is the boxlength in the i-th direction and the sum runs over the three- momenta q = 2k{j±, j^, £|) for 
integer rij. Following Ref. [20| we find 

i 

shAL) = E (w) 2 K i (V^) . ( B2a ) 

«yj ,i(£)=0, (B2b) 

1 roc . 

^U(£) = ^2 / dw E K ° ( a/I^R) e mp ' , (B2c) 

where l n = (riiLi, 77.2X2, n?,L-$) T , denote modified Bessel functions of the second kind and we have defined 

f 2 (u) ee uM 2 + (1 - u)m 2 + (u 2 - u)p 2 . (B3) 
For practical implementation we truncate the sum at n, = ±4 (i = 1, 2, 3). 
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Appendix C: Results by diagram 



In the following we present the results for the individual diagrams shown in Fig. Q] (at infinite volume) . For the 
matrix element of the Ioffe current they read: 



Amp* = \^Zia { °'N(k) , 



Amp[ = 16a¥ } mlN{k) , 



Amp). = — 



(0) 



3a\ 

(0) 



h,oN(k) 



i 3a\ 'g A 



[lip + {m 2 N - mg)/i : i - m N (m N + mo)4|ij N(k) 



(Cla) 
(Clb) 

(Clc) 
(Cld) 



(n) (n) 

One obtains the amplitudes for the Dosch current by the replacement a\ —> 2(i k . The amplitudes of the Chernyak- 
Zhitnitsky current read: 



Ampf = VZnf ^iz ■ k)N(k) , 



Amp^ z 



4/cf ] ml${z ■ k)N{k) 



AmpP 
Ampf 



19*<°> 
M, 



h fl t{z-k)N(k) , 



(l 1>0 + (m% - m 2 ) (/!,! - J«) - m w (miv + mo) (/g - ig)) ■ fc)2V(A) 
The Z-factor contribution has the form: 

= 1 + 7^4— (m^ + m N ) 7i,i + ^ I 



8F2 



(C2a) 
(C2b) 

(C2c) 
(C2d) 

(C3) 
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